CONTINUED FRACTIONS ON THE HEISENBERG GROUP 
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Abstract. We provide a generalization of continued fractions to the Heisen- 
berg group. We prove the ergodicity of the corresponding Gauss map and the 
convergence of the infinite continued fraction including an explicit estimate on 
the rate of convergence. Along the way, we show several surprising analogs of 
classical formulas about continued fractions. 



1. Introduction 

A regular continued fraction (RCF) expansion represents an irrational number x € 

E as 

1 

(1.1) x = ao H , ao € Z, <Zj € N,i > 1. 

ai 



a 2 H 

The integers CF{x) := {ao, ■ ■ ■ , } are the continued fraction digits of x (also called 
the partial quotients). Regular continued fractions and their many variations have 
played an important part in Diophantinc approximation, hyperbolic geometry, and 
the study of quadratic irrationals. 

Many higher-dimensional generalizations of continued fractions have been devel- 
oped to extend this powerful theory, but these efforts have met with varying success. 
In this paper, we develop a notion of continued fractions in the non-commutative 
setting of the Heisenberg group (in a sense, a complex two-dimensional contin- 
ued fraction). Surprisingly, we recover not only standard results of convergence 



and ergodicity (see Theorems 1.3 and 1.4), but also several simple, direct analogs 



of classical formulas for regular continued fractions — formulas which lack simple 
analogs for any other known multi-dimensional continued fraction. This suggests 
that continued fractions are a reasonable and natural object of study on the Heisen- 
berg group. 

This paper provides the basic properties of Heisenberg continued fractions, and 
opens up the way for many new questions. In future papers, we intend to link our 
study to that of Diophantine approximation on the Heisenberg group (see [7]) and 



the behavior of geodesies in complex hyperbolic space (see Remark 2.8 and |16|). 
Additional interesting questions include extending these results to similar spaces 
and their lattices (specifically, we expect our results to hold for all boundaries of 
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hyperbolic spaces) , a characterization of periodic continued fraction expansions, a 
careful analysis of the Gauss-Kuzmin problem, and a description of the dual space. 

The setting for this paper will be the Hciscnbcrg group EI, arguably the most natural 
non-commutative generalization of K. Specifically, H is M 3 with the modified group 
law (which we denote by *) 

(1.2) (x, y, t) * {x 1 , y', t') = (x + x , ,y + y', t + t' + 2(xy' - yx')). 

Note that in the first two coordinates one sees the usual addition of vectors, while 
the third coordinate incorporates an antisymmetric term. Note also that the group 
inverse (x, y, t) -1 of an element (x, y,t) € H is given by (—x, —y, —t). 

Let H(Z) denote the set of points in H with all integer coordinates. These form a 
subgroup of H, and we will think of them as the integers within H. Likewise, we 
think of points with all rational coordinates, H(Q), as rational points. 

Given a generic point ftgl there is a unique nearest Heisenberg integer [h] £ H(Z), 
with respect to the Heisenberg group's standard gauge metric: 

(1.3) \\(x,y,t)\\ = V(x 2 + V 2 ) 2 + t 2 d(h,k) = \\h-Uk\\. 

Note that left translations by elements of EI are isometries. That is, d(g *h,g*k) = 
d(h, k) for all g,h,k € EL In addition, one has an inversion operation (see j |2.1[ ) 
l : H\{0} -> H\{0} satisfying 

\\t{h)\\ = \\h\\-\ 

Given a point h £ H, we may remove the integer part of h via * h. We may 
then invert to uncover additional digits. 

Definition 1.1. The continued fraction digits CF(h) = {7^} and forward iterates 
{hi} of a point h £ EI are defined inductively by: 

70 = [h] h = 7Q -1 * h, 

Ti+i = Hh)] h+i = 7^\ * t(hi). 

Note that t(0) is undefined. Thus, the process may terminate after finitely many 
steps. We will characterize points for which this happens in Theorem |3.10| and, for 
the majority of the paper, focus our attention on points with infinitely many digits. 
We will also generally assume that 70 = unless otherwise specified. 

Definition 1.2. Let {ji} be a sequence of elements of EI(Z). For a finite sequence, 
define the associated continued fraction, 

K{7j} = K{7 4 }™ =0 := 70^71 1 • • • Lj n , 

supressing product notation and parentheses. It is clear that if CF(h) is finite, 
then KCF(h) = h. 

For an infinite sequence, we write 

K{ 7 i} = K{ 7 J- := lim K{ 7i }? = o, 

n—>oo 

provided the limit exists. 

Our main result is to show that K and CF define a valid notion of a continued 
fraction expansion for a point in EI. Namely, we prove 
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Theorem 1.3. Let h £ H and {7^} a sequence of elements o/H(Z). Then: 

(1) if HTjII > 3 for each i, then K{7i} exists and is unique regardless of whether 
{ji} is finite or infinite (Theorem 3.1). 

(2) A point h £ EI satisfies h — K{7i}™ =0 if and only if h £ H(Q) (Theorem 



3.10). 



(3) Every point in H /ias a continued fraction expansion. That is, for all h £ 



the limit WLCF(h) is unique and equal to h (Theorem 3.21). 



Throughout Sj3j we obtain variants of classical continued fraction results. We show 
a relationship between the denominator of a rational point and the length of its 
continued fraction expansion in Theorem 3.11| We find a recursive formula for 



the approximants K{7i}" =1 in Theorem 3.18 
h £ 



and show that the distance between 



and its approximants K{7i}™ =1 satisfies a variant of a classical relation in 
Theorem 3.23 We prove that the convergence of KCF(h) is uniform in Theorem 
1X261 

In fj4] we consider a generalization of the classical Gauss map x 1— > 1/x — |1 /x\ . 
Namely, let Krj C H be the Dirichlet region for H(Z), defined as the set of points 
h such that [h] — 0. It is easy to see that Kb is a fundamental region for H(Z), 
that is, the translates of Kq by elements of H(Z) tile H without overlap. It is also 
clear that for all h £ H, one has [h]~ 1 h £ Kb- 

We define a function T : Kb — > Kr> on the Dirichlet region by T(h) — [t/i] _1 i/i. 



We prove (see Theorem 4.101 



Theorem 1.4. There exists a measure fi on Kr> such that: 

(1) is T -invariant. 

(2) /j, is absolutely continuous with respect to Lebesgue measure on 

(3) T is ergodic with respect to /i. 



More generally, Theorem |4.10| implies the existence of a Gauss-Kuzmin measure 
on any fundamental domain K for H(Z) whose closure is strictly contained in the 
unit sphere (with appropriately modified Gauss map T). In particular, Theorem 
1.4 will provide a new measure on the unit cube Kc — [—1/2, 1/2) 3 , which is also 
a fundamental domain for H(Z). See Figure [4] for an illustration of the measures 
on Kc and Kp. 

We will now recall some background on classical continued fractions (j ]l.l[ ) and the 
Heisenberg group (J5| , and then study Heisenberg continued fractions in f|3] and 
the Gauss map in JlT 



1.1. Classical Theory of Continued Fractions. There are many variants on 
classical continued fractions and many ways to approach them (for good general 
references, see [51 O [HI HI])- We shall examine Nakada's a-continued fractions, 
since the study of them bears the most immediate resemblance to the Heisenberg 
continued fractions we examine in this paper. The a-continued fractions have 
two well-known continued fraction variants as special cases: Regular Continued 
Fractions (when a = 1) and Nearest Integer Continued Fractions (when a = 1/2). 
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T a x . — 



Let a € (0, 1]. Define the the a-Gauss map T a : [a — 1, a) — > [a — 1, a) by 

^0, x = 0, 

where [x] a is the unique integer such that x — [x] a G [a — I, a). Most continued 
fractions variants begin with these three simple pieces: a fundamental domain 
([a — I, a) here), an inversion that takes a point out of the fundamental domain 
(a; -1 ), and a linear translation that shifts us back into the fundamental domain 
Hz" 1 ]*). 

The digits of the a-continued fraction expansion for a number 16 [a — 1, a) consist 
of two parts, (a n ,e n ), where 

a n = a n (x) = [T^~ 1 x] a and e„ = e n (x) = sgn(7^ _1 a:). 

The sequence of digits (a„,e„) terminates when T™x = 0. These digits serve to 
record the data that is lost by iterating the non-injective map T a . In particular, 
we have 

ei ei 



x = 



ai 



a 2 + T^x 

Note that (a n {x), e n (x)) = (a n -i(T a x), e n -i(T a x)), so that T a acts as a forward 
shift of the continued fraction digits of x. 

One of the fundamental objects of study in the field of continued fractions is the 
sequence of convergents or rational approximants for a number x, given by 

Pn e _l 

In ' e 2 

ai H 



a 2 _| 1 

a n 

It is often easier to understand abstract properties of the sequence of convergents 
for a number x, than it is to understand abstract properties of the whole continued 
fraction expansion for x. 

A particularly useful property of convergents is the following matrix relation: 
Pn~l Pn \ _ ( ei \ f e 2 \ ( e, 



1 q n -i q n J \ 1 ai J \ 1 a 2 J \ 1 a„ 

From this relation, one can derive the recurrence relation q n = a n q n -\ + e n q n -2- 
While it would be nice to know that the q n form an increasing sequence of positive 
integers, this is not always the case (such as with continued fractions with odd 
partial quotients pQ). 

We can treat matrices as Mobius transforms, via 

a b \ _ az ■ 
c d 

If we do this, then the simple relation 
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together with (1.4 1, implies the more interesting relation 



(1.5) x = { Pn ^ Pn ) T£x = P^ T > + Pn 



q n -i Qn J q n -iT£x + q n 

By solving for T™x (or by applying the inverse of the matrix to both sides), one 
can obtain 

(1.6) Tlx = (-1) • qnX - pn , 

q n -i%-Pn-i 



Careful — but elementary — manipulation of the formulas (1.5| and (1.6 1 yields 

n 

(1-7) q n x-Pn = (-l) n U T > = (-1)" ' ~ ^+i'" £ " 



i=0 



Tin -hi i 
a x + q n+1 



From ( 1.7| it is short exercise to see that q n x — p n converges to 0, and hence that 
Pn/q-n converges to x. Thus it makes sense to write x as an infinite continued 
fraction expansion 



£2 

Ol 



a 2 



There are varying notions of convergence for continued fractions variants besides 
the fact that \x— {p n /qn)\ tends to 0, which is typically known as weak convergence. 
In multi-dimensional continued fractions, where one might have convergents 



Pl,n P2,n Pk,n 

q n ' q n q n 



to a point (xi, X2, ■ ■ ■ , Xk) , 



the property that \q n Xi — Pi lU \ tends to for all i is known as strong convergence. 
(The Jacobi-Perron continued fraction, which is in many ways considered to be the 
prototypical multi-dimensional continued fraction, does not satisfy strong conver- 



gence.) The fact that all columns of the matrices ( 1.4 1 converge (projectively) to the 
same point is known as uniform convergence. Uniform convergence is non-trivial 
for higher-dimensional continued fraction variants. 

In general, it is hard to know whether an arbitrary sequence of continued fraction 
digits (a„,e„) £ R 2 produces a convergent infinite continued fraction. (Even the 
seemingly innocuous two-digit sequence {(1,1), (1,-1)} causes convergence prob- 
lems.) One major result on this question is Pringsheim's theorem, which states 
that if \a n \ > |e n | + 1 for all n £ N, then the infinite continued fraction converges. 
For more on this topic, see [TT] . 

For many continued fractions, the digit shift map T is ergodic with respect to some 
invariant measure. For Regular Continued Fractions (a = 1), the invariant measure 
that is absolutely continuous with respect to Lebesgue is the classic Gauss measure 

i r i 

u(A) = - — - / dx. 

^ ' \og2j A l + x 

The ergodicity of the map T means that there is a notion of average behavior for 
continued fractions. For example, for almost all x, the regular continued fraction 
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expansion of x satisfies 

logg„(x) 7T 2 



lim 



n 12 log 2 

and a n — 1 approximately 42 percent of the time. 

Applications of continued fractions come from various areas. We mention only a 
few in greater detail here. One of the most classical results on continued fractions 
is Lagrange's Theorem, which states that x has an eventually periodic continued 
fraction expansion if and only if a; is a quadratic irrational number: thus one often 
studies properties of quadratic irrationals by understanding their RCF expansion. 



The term q n x — p n that appeared in (1.7) is closely related to the study of best 



approximants — namely, rational numbers n/m that satsify the following relation 

\mx — n\ < \bx — a\ a,b £ Z, 1 < b < m 
must be an RCF convergent p n /q n for a;. 



2. The Heisenberg Group 



We will think of the Heisenberg group in three different ways. For geometric pur- 
poses, including illustration and discussion of measures, we will identify EI with 
M. (with the appropriate group structure and geometry). For the majority of the 
paper, however, we will be concerned with the representation of H as a group of 
unitary matrices or as a subset of C 2 . This is in direct analogy with thinking of 
the real numbers as elements of SL(2,M.) or as the real axis within C 1 . We now 
discuss these models, and then record some information on discrete subgroups of 
H and their fundamental domains. 

We emphasize that the topological and measure-theoretic notions we consider do 
not (qualitatively) depend on the model we choose, nor on the metric. In particular, 
convergence in H can be shown using the intrinsic gauge metric, or using metrics 
intrinsic to the model, such as the Euclidean metrics on R 3 or C 2 . 

2.1. Geometric Model. In the introduction, we defined H as the space R 3 with 
group law 

(x, y, t) * {x 1 , y', t') = {x + x',y + y', t + t' + 2(xy' - yx')). 

Combining the first two coordinates into a complex number, H becomes CxR with 
group law 

(z, t) * (z', t') = (z + z', t + 1' + 2 Im(zz')). 

We will think of these as the same model, and use it primarily when geometry or 
visualization are concerned. There are several standard (topologically equivalent) 
metrics on H; we will work with the gauge metric. The gauge ||-|| and distance d 
are defined by: 

||(M)|| = \f\z\ 4 + tz d(h,k) = \\h- 1 *k\\. 

There are four basic transformations we are interested in: 
(1) Left translations h h-> k * h, for k € H, 
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Figure 1. Two views of nested spheres in H, centered at (i,j,0) 
with i,j 6 {— 1,0,1}, related to each other by left translation by 
elements of H(Z). In the top view (left), the spheres look identical. 
A side view (right) shows an additional a shear in the t coordinate. 



(2) Rotations z n- e i6 z, for eR, 

(3) Metric dilations S r : (z,t) M> (rz,r 2 i), for rel, 

(4) The Koranyi inversion i : H\{0} — > H\{0} given by 



i(z,t) 



it 



t 2 



Translations and rotations do not distort distances or volume (that is, the Lebesgue 
measure A on IR 3 ). The map S r is a group homomorphism dilating distances by a 
factor of r and volume by a factor of r 4 . The Koranyi inversion is a conformal map 
with the following important property. 

Lemma 2.1 (See p.19 of 0). Let h,ke H\{0}. One has 

d(h,k) 



d(ih, ik) 



\h\\\\k\\ 



In particular, one has \\ih\\ = \\h\\ 1 , so that the inside and outside of the unit ball 
are interchanged. Note that individual points on the unit sphere are not fixed. 



Remark 2.2. We will show in Lemma 2.12 that i has a particularly simple form 



in the unitary model. It is also conformal with respect to the gauge metric [5] . 

We record the following relationship between volumes and radii of balls in H. In 
particular, the lemma implies that the Heisenberg group has Hausdorff dimension 
4, and that A is equivalent to the Huasdorff 4-measure on H. 

Lemma 2.3. The volume of a ball B(h, r) of radius r around a point h is given by 

X(B(h,r)) =r 4 A(B(0,l)) 
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Proof. Applying a left translation, we may assume h is the origin. Further, we may 
rescale 5(0, r) by the Heisenberg dilation 5u r (z,t) — (z/r,t/r 2 ) to obtain B(0, 1). 
The dilation distorts A by a factor of r~ 4 . □ 




Figure 2. Spheres in H centered at the origin, with radius 
2,1,1/2, with sectors removed to display nested spheres. The 
spheres are parametrized by applying t to a plane; the radial lines 
of the plane provide the characteristic foliation on the spheres. 



2.2. Real Nilpotent Model. It is common to describe HI as the group of nilpotent 
upper-triangular 3-by-3 real matrices. Our definition is related to this real nilpotent 
model via the Lie group isomorphism 




(2.1) {x,y,t)^ 
We will not use the real nilpotent model, although our results can be rephrased for 



it. Note that under (2.1 ), H(Z) is not identified with matrices with integer entries. 



2.3. Unitary Representation. For calculation purposes, we will use the (Siegel) 
unitary representation of EL Namely, we will embed EI in GL(3,C) via the homo- 
morphism: 

/ 1 

(2.2) U:(*,t)i-> z(l+i) 1 

\ \z\ 2 +ti z(l-i) 1 

Remark 2.4. In literature, one sees a factor of v2 rather than 1 + i in the em- 
bedding. The latter is more convenient for our purposes. 
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Let J be the Hermitian inner product given by 

H( z o,Zi,z 2 ),(w ,w 1 ,w 2 )) = ( 2o zi Z2 ) 
In particular, we record 

(2.3) \(z ,zi, z 2 )|j = S((z ,z 1 ,z 2 ),(zo,z 1 ,z 2 )) = |zi| 2 - 2Re(z^z 2 ). 

We will refer to a vector of norm as a null vector. 

Abusing notation, we will also use J to denote the skew-diagonal matrix above. 
Note that J has signature (2, 1): it has two positive and one negative eigenvalue. 

The unitary group SU{2, 1) C GL(3, C) is the set of matrices M e SL(3, C) satis- 
fying J(Mz, Mw) = J(z, w) for all z, w E C 3 . Equivalently, M satisfies Aft JM = J, 
where f denotes the conjugate transpose. It is easy to see that: 

Lemma 2.5. U(H) C SU(2, 1). 

Later calculations will require us to step outside of U(H). The following lemma 
provides a basic property of elements of SU(2, 1). 

Lemma 2.6. Every matrix in SU(2, 1) is of the form 




02,3 a l,l — 02,101,3 
(2.4) ( a 2 ,i 03,301,1 ~ Q3,iQi,3 

03,302,1 — 03 ; i(Z2,3 

Proof. Every matrix M = {a id ) in SU(2, 1) satisfies M f J = JM -1 . Wc have 




Aftj 




02.1 -Ol,l 

02.2 — «1,2 

02.3 — a l,3 



On the other hand, 



a 3,io 2 2 — 03,202,1 03,2^1, 1 — 03 ; iai ; 2 02,101^2 — 02,201,1 
JM -1 = ( a 3i ia 2 ,3 - 03,302,1 03,301,1-03,101,3 02,101,3-02,301,1 

03,202,3 — 03,202,2 03,301,2 — 03,201,3 02,201,3 — 02,301,2 

Comparing the two matrices completes the lemma. □ 

2.4. Siegel Model. The Siegel model provides a geometric view of the unitary 
representation and a simpler formula for the Koranyi inversion. We will in fact 
define two closely related models, the planar Siegel model that views a point h&M 
as a vector (u, v) G C 2 , and the projective Siegel model that views h as a point in 
complex projective space with homogeneous coordinates (1 : u : v). We will denote 
both models by S. 

We first identify a point h £ H with geometric coordinates (z,t) with the vector 

(2.5) (l,z(l + i),|z| 2 +it) eC 3 . 

Note that this is exactly the image of the vector (1,0,0) under the unitary trans- 
formation U(z,t). We will say that h has planar Siegel coordinates 

(2.6) (z(l+i),|z| 2 +it) e C 2 . 
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The planar Siegel model of H is the set of points in C 2 of the form (2.6 1 



Some times, a unitary transformation will take (1, z(l \z\ + it) to a point that 
is not of the same form, but can be rescaled to be such. It will therefore be useful 
to think of vectors up to rescaling, that is, as elements of complex projective space 
CP 2 . 

Recall that the complex projective plane CP 2 is the projectivization of C 3 , i.e. the 
set of non-zero vectors up to rescaling by a non-zero complex number. A point in 
CP 2 has homogeneous coordinates (zq : z\ : Z2), well-defined up to rescaling. 

We can now define the projective Siegel model of H as the set of points in CP 2 with 
homogeneous coordinates (1 : z(l + i) : \z\ 2 +it). 

Abusing notation, we will denote both Siegel models by S, with the identification 
(it,?;) ■<->• (1 : u : v). We have the following simple characterization of points in S. 

Lemma 2.7. Let (zq : z\ : z 2 ) £ CP 2 be a null point, that is \\(zo, Z\, ^2)||j = 0. 
Then either (zq : z\ : z 2 ) € S or (zq : z\ : z 2 ) = (0:0:1). 

Remark 2.8. The region {(z : z\ : z 2 ) € CP 2 : || (z , zi, z 2 )||j < 0} bounded by 
S U {(0 : : 1)} is the Siegel domain. Complex hyperbolic space is defined on this 
region and has strong connections to the Heisenberg group, see e.g. [21 03 [HI 02] ■ In 
particular, we intend to discuss the relation of Heisenberg continued fractions to 
geodesic coding in complex hyperbolic space in an upcoming paper, following [15] . 

Note that the gauge norm is easy to write in the Siegel model: 

1/2 

Lemma 2.9. Let (u,v) € S. Then the gauge norm of (u,v) is \\(u,v)\\ — \v\ . 
Proof. An clement of S has the form {u, v) = (z(l + i), |z| 2 + ti ) for some (z, t) £ 



H. The gauge norm of (z,t) is given by ||(z,t)|| = d\z\ A + t 2 — \v^ 2 . □ 

The gauge distance is defined as d{h, k) = ||/i _1 A; . With this in mind, we show: 
Lemma 2.10. In the planar Siegel model, we have 

{ui,v 1 )~ 1 {u 2 ,v 2 ) = (u 2 - Ux,W[-uiu 2 +v 2 ). 

Proof. We have associated to (ui,wi) _1 the matrix 

1 

-Ui 1 

«1 —Ml 1 

Applying this matrix to the point (1,1*2, 1*2), we get the vector 

(l,1t 2 - Ux,V\-UlU2 + V 2 ). 

Taking the last two coordinates yields the desired formula. □ 

We now study the action of SU(2, 1) matrices on the Heisenberg group in the 
Siegel models. General linear matrices act on CP by acting on the homogeneous 
coordinates. Since we have C 2 CP 2 by taking (u, v) h-> (1 : u : v), we also obtain 
an action on C 2 . 
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Lemma 2.11. Let M = (cnj) e Gi(3,C) and (u,v) eC 2 4 CP 2 . Then M acts 
on (u, v) as: 



M(u,v) = 



02,1 + ^2,2" + &2,3V 03,1 + a 3,2" + 03,3V 
ffll.l + 0>1,2U + Oi,3« ' 01,1 + Oi^M + 01,3V 



Proof. The point (u, v) corresponds to a point in CP 2 with homogeneous coordi- 
nates (1 : u : v). We then have 

ai,i + ai i2 w + ai,3« 

02,1 + 02,2« + a 2 ^v 

03,1 + a 3 ^u + a 3t3 v 

To view M(l : 11 : n) as a point in C 2 , we renormalize so that the first coordinate 
is 1, and take the remaining two coordinates. □ 





( 1 > 




M 


u 


-( 




{ v J 





Elements of GL(3, C) do not necessarily preserve the set S, but the unitary matrices 
SU(2, 1) preserve J and therefore S. In particular, elements of U(H) act transitively 

0-1 

on S while fixing the point (0:0:1). Denote the matrix I 1 I by 



-10 



U( t ). 

Lemma 2.12. U(t) acts on H 6y i/ie Koranyi inversion t. 



Proof. We compute, for a point in H with geometric coordinates (z, t) and projective 
Siegel coordinates (1 : z(l + i) : \z\ 2 + ti): 

U(i)(l : z(l + i) : \z\ 2 + ti) = (\z\ 2 + ti : -z{\ + i) : 1) 

-(u 



= 1: 



|z| + ti 

— z 



ti 



T—ri 1 + 1 ) : m 



ti 



+ t 2 



= 1: 



U +ti 



-(1 + i): 



\zf+ti 



+ 



z\ 4 +t? 



We thus have that under U(t), the geometric coordinates are mapped to 

(w^+ii'ra)' asdesired - D 



2.5. Lattices and Fundamental Domains. Recall that H(Z) is the set of Heisen- 
berg points with integer coordinates. In the geometric model H = C x I, we 
have H(Z) = Z[i] x Z. In the Siegel model, H(Z) is the set of points (u,v) such 
that u, v <G Z[i], and w has norm divisible by 2. In the unitary model, we have 
H(Z) C SU(2, 1; Z[i]), where the latter denotes the set of unitary matrices with 
Gaussian integer coefficients, and is known as the Picard modular group. 

Likewise, we will denote by H(Q) the set of points in H with rational coordinates. 
Recall that the Hciscnberg group admits a family of dilation maps 6 r given by 
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5 r (z,t) = (rz,r 2 t) in the geometric model. The dilation maps are group isomor- 
phisms and satisfy d(S r h,S r q) = r ■ d(h,k) for all h, k G H and r > 0. It is clear 
that h G H(Q) if and only if there is an integer n€N such that S n h G H(Z). 

We are now interested in the structure and geometry of H(Z). We record its gen- 
erators in the geometric model: 

Lemma 2.13. The group H(Z) is generated by the elements (1,0), (i,0), and (0,1). 

As Falbel-Francics-Lax-Parker showed, Z[i] and SU(2, l;Z[i]) are closely linked: 

Theorem 2.14 ([!]). The group SU(2, l;Z[i]) is generated by the matrices U(1,0), 
U(0, 1), V(l), and il , where I is the identity matrix (see I 2.3\2.4 for the notation). 



Remark 2.15. Note that we defined SU(2, 1; Z[i]) with a particular Hermitian 
form J in mind. Different Hermitian forms J provide isomorphic Lie groups U(2, 1), 
but the lattice SU(2, l;Z[i]) depends on the choice of the Hermitian form. If two 
forms are related by an integer change of coordinates, then the associated lattices 
are equivalent. If the change of coordinates is not integral, the lattices are not 
isomorphic as groups (even up to finite index) , see [131 114) Nonetheless, in literature 
one mostly sees mention of the Picard modular group, defined by a Hermitian form 
equivalent to our J. 

We now discuss fundamental domains for H(Z). Recall that a fundamental domain 
for H(Z) is a connected set K c H with piecewise smooth boundary whose trans- 
lates tile H without overlap. That is, \j{g*K : g G H(Z)} = H and Kf] {g*K) ^ 
implies g = 0. 

Our definition of fundamental domain is slightly stronger than usual. The standard 

o o o 

definition only considers K n (g * K) ^ 0, where K is the interior of K. We require 
"half-closed" fundamental domains, so that the following holds: 

Lemma 2.16. Let K be a fundamental domain for H(Z). Then the map [p] K 
mapping all points of gK to g is well-defined. 

The following lemma follows immediately from the definitions: 

Lemma 2.17. The following regions are fundamental domains for H(Z): 

• The unit cube K c = [-1/2, 1/2) x [-1/2, 1/2) x [-1/2, 1/2). 

• The Dirichlet domain Kr> = {h G H : d(0, h) < d(g, h) for all g G H(Z)}, 
with a choice of excluded boundary points. 

Denote the unit sphere in H by S. For a subset A C H, let rad(^4) denote the 
supremum of the norms of the points of A, and let X(A) denote its Lebesgue measure 
(in the geometric model). 

Lemma 2.18. Every fundamental domain K for H(Z) satisfies X(K) = 1. Fur- 
thermore, the domains Kq and Kr> satisfy rad(Xc) = rad(i^D) = \/l/2. 

Proof. The radius of Kc is easy to compute because ||-|| behaves similarly to the 
Euclidean norm. As in the Euclidean case, the norm is maximized by each corner 
of the cube. We have |j (1/2 + il/2, 1/2)|| = {flj2. 



HEISENBERG CONTINUED FRACTIONS 



13 




The radius of Kd seems difficult to compute directly, as the boundary of Kd is 
more complicated (see Figure [3]). We will therefore argue indirectly by means of 
Kc- Let h £ Kd, and choose g £ H(Z) so that g * h £ Kc- We then have 
\\g*h\\ < ra,d(K c ) = This implies that d{g~ l ,h) < 0J2. Now, by 

definition of Kd, d(0,h) < d(g^ 1 ,h) < y/l/2, so we must also have \\h\\ < \/l/2. 
To prove equality, one shows that the point (1/2 + if/2, 1/2) is contained in Kd- 

For the volume computation, it is clear that X(Kq) — 1. To compute X(K) for 
an arbitrary fundamental domain K, note that Lebesgue measure is preserved by 
left translation in the Heisenberg group (which acts by shears). Since Kc can be 
constructed by rearranging measurable pieces of K, the two fundamental domains 
must have the same volume. □ 



3. Heisenberg Continued Fractions 
Fix a fundamental domain K for the group H(Z) such that rad(i^T) < 1 (e.g. K, 



c 



or Kd in Lemma 2.17). We begin by establishing some notation. 



Definition 3.1. Given an arbitrary sequence {7i}?=i 01 non-zero digits in H(Z), 
we write the associated continued fraction as 

(3.1) K{ 7l } = K{ 7i }?=i = niH2 ■ ■ ■ Hn- 

For an infinite sequence {ji}^, we define = K{^ i }'*L 1 := lim^oo K{7i}" =1 , 

if this limit exists. The goal of this section is to show that the limit does exist in 
several important cases, and that the computation of K{7i} may be simplified by 
using a recursive algorithm. 



Definition 3.2. We associate with K: 
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(1) A "nearest-integer" map [■] : H — > H(Z), characterized by 

[ft] = g for each g £ H(Z) and ft £ gK. 

Note that [•] selects the nearest integer in the gauge metric exactly if K is 
the Dirichlet domain Kjj. 

(2) The Gauss map T : K\{0} K given by 

T(ft) = [Lh}-hh. 

Remark 3.3. Working with the geometric model, one sees that each axis is pre- 
served by the Gauss map T. In particular, the restriction of T to each axis is 
essentially isomorphic to the nearest-integer Gauss map on [—1/2, 1/2]. The theory 
of continued fractions we develop likewise restricts to the classical nearest-integer 
continued fraction theory on the axes. 

Definition 3.4. Given a point ft £ K, have: 

(1) The forward iterates hi := T l h £ K, for each i, 

(2) The continued fraction digits 7.; := [ift^-i] £ H(Z), for each z, 

(3) The rational approximants K{7i}" =1 £ H(Q), for each n. 

Because T is defined on if \{0}, the process of defining forward iterates, continued 
fraction digits, and rational approximants terminates if for some i we have hi = 0. 
We will characterize the points h for which this happens in Theorem |3. 10 

More generally, for a point fcelwe can take 70 = [ft], fto = [7cT 1 ft] an d obtain the 
remaining digits {ji}°^ 1 of CF(h) from h$ £ K as before. However, our focus will 
be on points in K. 

It is easy to see that, on finite sequences, K is the inverse operation to CF: 

Lemma 3.5. For h £ K with CF(h) a finite sequence, we have KC-F(ft) = h. 

Remark 3.6. The operation K is defined without reference to a specific fundamen- 
tal domain K. Thus, while we will show that K.CF(h) — ft, we do not in general 
have CF(K{^fi}) = {7^}. Indeed, problems arise when the ji get too close to the 
unit sphere. 

For example, let K = Kc, the unit cube, and {7^} = {(01, &i) = (1, 0)}. We have 

K{ 7i } = t(l ) 0) = (-1,0). 

Attempting to reverse the process, we have (a ,b ) = [(—1,0)] = (—1,0), and 
(-1, O)" 1 * (-1, 0) = (0, 0), so that 

CJ?(-1,0) = {(00,60) = (-1,0)}. 

This non-uniqueness of continued fraction expansions is analogous to how in regular 
continued fractions we have, for example, 

1 1 
5 + i ~ 6' 
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3.1. Pringsheim-Type Theorem. The Pringsheim Theorem for regular contin- 
ued fractions guarantees convergence of a continued fraction whose digits are suffi- 
ciently large. A variant holds for the Heisenberg group: 

Theorem 3.7 (Pringsheim-Type Theorem). Let be a sequence of points 

in H(Z) such that for each i we have \\"fi\\ > 3. Then the limit K{7,;} exists. 
Furthermore, CF(K{7;}) = {7^}. 



Proof. Recall that left multiplication by an y 7 <£ H(Z) is an isometry, and l satisfies, 
for all h,k e 



d(ih, ik) = (Lemma 



2.121 



Let Ko be the Dirichlct fundamental domain for H(Z). It follows from the definition 
of Kd and the triangle inequality that for each point h <E H with \\h\\ < 1/2, we 
have h £ Ku- Conversely, for each point h € Kr> we have by Lemma [2.18| that 

IWI<</T72. 

Suppose that 7 G H(Z) with ||7|| > 3. We claim that i^Kd C Kd- Indeed, every 
point h € jKjj satisfies \\^fh\\ > 3 — \\h\\ > 3 — \A/2 > 2, so that \\i/~fh\\ < h, and 
we conclude l'jKd c Kjj. 

Now, for each n, we have (because the identity element is contained in Kjj): 

K{7i}™ =1 = 471672 • • • ijn 
= niil2 ■ ■ ■ ^7n0 
<E 471.172 • • • HnK D . 

These cylinder sets form a nested sequence: 

LllH2 ■ ■ ■ HnK D = t7lt72 • • • Hn-l(HnK D ) 

c 171172 • • • fy n -iK D - 



By the above calculation, the diameter of the cylinder set 471^72 ••• Lj n Krj is 
bounded above by (3 — ^/l/2)^ 2n diam(Kr)). We thus have that the sequence 
of fractions K{7i}™ =1 (as n varies) is a Cauchy sequence, and hence converges to 
some K{7i}. 

We thus have that exists. By construction, we also know that it is contained 

in the cylinder sets ljil^ ■ 1 1 l 1uKd, for each n (note that the cylinder sets are in 
fact properly nested, so that cannot escape to a cylinder set's boundary). 

This is equivalent to the second assertion of the theorem. □ 



3.2. Rational Points. We will now show that a point in H has rational coordinates 
if and only if it has a finite continued fraction expansion. Our proof is motivated 
by the work of Falbel-Francsics-Lax-Parker [3] and uses the Siegel model. 

Recall that for a point h € K that is of interest to us, we will write 

h = (u,v) £ C 2 

in the planar Siegel model. We will also think of (it, v) as the element of CP 2 with 
homogeneous coordinates (1 : u : v). That is, it is the vector (l,u,v) considered up 
to multiplication by a non-zero compex number. 
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Definition 3.8. Given an element 7 <E H(Z) with planar Siegel coordinates (a, /?) <E 
(Z[i] xZ[i])n5, define 

/ 0-l\/l0 0\ 
A 7 := U(t)U(7) = 1 a 1 

\ -1 / \ p a 1 / 




Lemma 3.9. In £/ie Siegel projective model, we have 

K{ 7i }?=i= A^.-A^ (1:0:0). 



Proof. Abstractly, we have the definition K{7i}™ =1 = 471 1 • • • i7„. Using the identity 
element £ H, we may also write K{7i}™ =1 = L-fiL ■ ■ ■ L-f n 0. In the projective Siegel 
model, is interpreted as the point (1 : : 0) £ CP . The inversion 1 and left 
multiplication by 7$ are, respectively, interpreted as the unitary matrices U(t) and 
U( 7i ). Thus, K{ 7i }? =1 = A 7l • • • Ay n {l : : 0), as desired. □ 

We are now in position to characterize rational Heisenberg points in terms of their 
continued fraction expansion. 

Theorem 3.10. Let h £ H. Then h £ H(Q) if and only if h = K{ 7l }f =0 for some 
finite sequence {"fi}" =0 - 



Proof. Suppose h = K{7i}" =0 
fact that 7 l £ H(Z) that h £ 



Then it is clear from the definitions of K and the 



Conversely, fix K = Kjj and assume by way of contradiction that there exists an 
element h £ H(Q) with an infinite continued fraction sequence CF(h) = {7^}°^. 
Without loss of generality, we may assume h £ K (this corresponds to discarding 
the digit 70 of h). 

The idea of the proof is to show that the forward iterates hi of h can be written as 
fractions whose denominators decrease with i. Write, in planar Siegel coordinates, 

r p x 



with q, r,p £ Z[i], Because h £ K,we have by Lemma 2.9 that \p/q\ < r&d(K) 2 < 1. 
Consider the first forward iterate hi = Th = 7 1 _1 t/i as a vector in C 3 : 



(i) 



P 












(1 





1 




r 


-1 


-ST 




\ P 



-p 

r + aip 
-q - air - /3ip 



Thus, hi is a rational point with planar Siegel coordinates hi 
thermore, we have q^ = p, so that 



\qM ' g(l) 



Fur- 



(3.2) 







p 


q 




q 



\h\\ 2 < rad(if) < 1. 
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Repeating this procedure recursively, we have rational coordinates hi = f^jy, W>) 

for each forward iterate h i: satisfying — p^ 1 ^. Since h t G K for each i, we 
obtain for each n: 

(3.3) 



>) 



< M (rad(tf)) 



2n 



For sufficiently large n, we conclude |g„| < 1, which implies that q n 
is only possible if h n _ 1 = and CF(h) is, in fact, finite. 



0, but that 
□ 



As a corollary to the proof of Theorem 3.10 we obtain 

Theorem 3.11 (Denominator Growth Theorem). Let h G H(Q), wzift CF(/i) = 
{7i}™=0- Suppose one can write h as a fraction with denominator q G Z[i]. Then. 

\q\> 2"/ 2 . 



Proof. Th e resu lt follows directly from (|3.3|), using either fundamental domain K 
in Lemma 



2.18 



with radius bounded by \/l/2, and the fact that q n ^ 0. 



□ 



Remark 3.12. One may hope for a stronger statement that for a sequence {7i}^i 
of elements of H(Z), the norms of the denominators q n of the partial fractions 
K{7i}™ =1 are an increasing sequence. However, we are unable to prove this without 
assuming that ||7j|| > 2 for all i. 

3.3. Recursive Formula. We will now find a simple recursive formula for K{7i}^ 1 . 
Definition 3.13. Let {7^} be a sequence of elements of H[Z], Define 

Qn • — ^71 ' ' ' ) 

(q n ,r n ,p n ) := Q„(1,0,0). 
We have the following by Lemma |3~9] 

Lemma 3.14. In the above notation, K{7i}" =1 = f^N > ^ n ^ e planar Siegel 
model. 

Remark 3.15. It should be noted that Theorem |3. 1 1 does not imply that q n > 
2™' 2 . Recall from Remark 3.6 that if CF(h) — {~/i}f^L , we do not necessarily have 
that CF(K{ 7i }f =0 ) = { 7i KL . 



Lemma 3.14 states that the partial fraction K{7i}™ =1 is encoded in the matrix Q r 
As in the case of regular continued fractions, Q n stores additional information: 

Lemma 3.16. In the above notation, the matrices Q n have the form 

Qn — 



Qn 


q„ 


-Qn-1 


r n 


fn 


-r n -i 


Pn 


Pn 


-Pn-1 



where the elements q n ,f n ,p n are given by 

q, 

P 



r n q n - 


1 


qn^n— 1 1 


p n q n - 


1 


- q n Pn-l 


Pn^n- 


-1 


- r n p„_i 



18 



A. LUKYANENKO AND J. VANDEHEY 



Proof. The first column of Q n is as stated by the definition of the vector (q n , r n ,p n ). 
The third column follows from the identity Q n = Q n -\A ln . The second column 
follows from the "cross product" Lemma [2. 6| □ 



We record the following for later use: 

Lemma 3.17. The following formula holds for Q n : 



-Pn 


r„ 












-pn 


f n 


-Qn 










-Pn-1 


r n -i 


-Qn-1 J 












Priori 


Priori 


p n q n 


- q n p n 


™nQn 


' n Qn 


■( 


Pn-%r n 


-p n r n -i 


p n q n -i 


-Pn-XQn 


r n -iq n 


- r n q n -i 




Pn-lfn 


-p n r n -i 


p n q n -i 


- Pn-l<jn 


r n -iq n 


- f n q n -\ 



Proof. This follows from the identity QLJ = IQn 1 and Lemma 3.16 



□ 



We can now obtain a recursive form for the partial fractions K{7j}™_ 1 . 

Theorem 3.18. Let {7^}.^! be a sequence of elements o/H(Z) represented in the 
planar Siegel model by the vectors {{on, pi)}^.^. Let (q_i,p_i,r_i) = (0,0,1) and 
(qo,POi r o) — (1)0,0). Define, recursively, for n > 0, 



Qn+l ) 




( Qn 


r n +i 


- 


r n 


Pn+l ) 




\ Pn 



q n r n q n -i - q n r n -i -q n -i 



r„ p n q n ~i - q n Pn-i -r n -i 



Pn P n r n -1 ~ r n p n -l -Pn-l 



Then for each n we have, in the planar Siegel model, 




K{ 7l }?=i = 



T71 Pn 
qn Qn 



Proof. Earlier in the section, we defined matrices A li (which append the digit 7, to 
a continued fraction) and Q n = A yi ■ ■ ■ A ln . We set (q ni r n ,p n ) = Q n { 1,0,0). We 
claim that this agrees with the definition in the statement of the theorem. Lemma 

3.14|will then tell us that K{ 7 J™ =1 = 
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Taking Qq to be the identity matrix, the following computation provides the equiva- 
lence (see the definition of A^^, and Lemma|3.16|for the form of the two matrices) 





— Pn+1 


-Oin+l 






Oin+l 


1 


: 





-1 







\ o 



Qn 


q n 


-q n - 


1 


r n 


Pn 


-r n - 


1 


Pn 


r n 


-Pn- 


1 




Rewriting q n , f n ,p n in terms of the other terms in Q n completes the proof. 



□ 



3.4. Continued Fraction Representation Theorem. We are now ready to 
prove the convergence of continued fraction expansions. In fact, we obtain a varia- 
tion on the strong convergence property. While we do not obtain strong convergence 
in the sense of Schweiger |15j . our convergence estimate is obtained via a similar 
method to strong convergence for regular continued fractions. We hope to im- 
prove this estimate and explore applications to Diophantine approximation in an 
upcoming paper. 

We also note that we obtain such an explicit convergence estimate by exploiting a 
special form for Q^ 1 that follows from the identity M^JM = J that defines U(2, 1). 
Other continued fraction theories are complicated by the lack of a simple form for 

Before we can prove convergence, we need to show that q n will never equal 0. We 
prove this in two steps. 

Lemma 3.19. We have 



(3.4) 



q n + q n u n 

Pn + PnU n 




(-1)" 



V 



1 


VV\ 


■■Vn-1 




u 




••1>r»-l 




1 



Proof. By Lemma 3.16 the vector on the left-hand side of (3.4) equals 

Qn(lj u nj V n ) — Aj ± ■ ■ ■ A^ n (1, 1i n , V n ). 

Recall that the forward iterates of h are given by hi = T l h = A^ 1 ■ ■ ■ A~^h, and 
have planar Siegel coordinates (ui,Vi), corresponding to the points (1 : Ui : Vi) € 
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More generally, we have A li ■ ■ ■ A ln (\ : u n : v n ) = hi. Write A ln (1, u n , v„) 
(a, b. c). Since A~ n has the form (see Definition 3.- 



—a 




1 


:)■ 









-l,V n -l), 




U n -1 




l' «n-l 



we have that c = —1. Since {b/ a, cj 

^7,i ( 1 1 "in u n ) 

Continuing in the same fashion we see that 

Aj n _ i ^7„ (lj^iufre) = ^-7„ _ i ( j j — 1 

\ ^n— 1 Un— 1 

1 1 



After n iterations, the process yields the desired formula. □ 
Lemma 3.20. For n > 0, we have that q n never equals 0. 

Proof. Assume, by way of contradiction, that q n — 0. Then by Lemmas |3.16| and 



3.17 we have that q n = as well (r n also must equal 0, but we will not use this 



fact). Since the matrix Q n has determinant (—1)™ and each entry is a Gaussian 
integer, we must have — q n -i must have norm 1. 

Therefore, we have that 

(3.5) \q n + q n u n - q n -iv n \ = \v n \ < 1 
However by Lemma |3.19| we have that 

(3.6) \q n + q n u n - q n -\V n \ = \vV\V2 ■ . .w„-i| _1 > 1, 

which is a contradiction. Therefore our assumption that q n — must be false. □ 

Now we can continue with the proof of convergence. 

Theorem 3.21. Let h G H and let K be a fundamental domain for H(Z) with 
rad(JT) < 1. Then 

KCF(h) = h. 

Furthermore, if CF(h) — {7^} is a sequence with at least n terms, then the rational 
approximants satisfy 

d(K{ ll }? =0 ,h)<Y a d(K) n+1 

for both rational and irrational points in H. Let q n be the denominator of the n th 
rational approximate. Then we in fact have 



d(K{ 7i }? =0 ,/i)< 



1 1 1/2 
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Proof. Recall from Lemma 3.14 that the associated rational approximates K{7i}™ =1 
have planar Sicgel coordinates (jj^, <p)> associated to the vector (q n ,r n ,p n ) € C 3 . 
Recall also that the forward iterates T n h have planar Siegel coordinates (u n ,v n ), 

1/2 

and we have \v n \ < rad(if) < 1. 
It therefore suffices to show that 



d(KH}f =0 ,h) 



n 



i=0 



,1/2 



,1/2 



Indeed, by Lemmas 2.10 and 2.9 we have 



d(K{ 7i }? =0 ,h) = d( ( 



U , V ■ 

q n 



v — 



u - 



In 



1/2 



1/2 



We now view h as the vector (1, u, v) and represent the operation T n by the unitary 
matrix Q^ 1 - The vector 



Qn 1 



is then a scalar multiple of (l,u n ,v n ). In particular, 



1 




^ -Pn-l ?*«-l 


-q n -i 




I- 




-q n 






V Pn -r n 


q„ 







) 









Pn-i ~ r n ~iu + q n -iv 
By multiplying this formula together for various indices we obtain 

rt. n. 

Pi - TiU + qtv 



n«*=(- i ) n n 



Pi-i - n-iu + qi-iv 



8=1 

pT) - rf)U + cfiv 
■. n P^-r^u + q^v 



= (-1)' 

= (-ry 

V 

This yields the interesting formula 

71 

(3.7) p^-f^u + -q-v = (-l) n Y['< 

i=0 
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Wc then have 



d(K{ 7i }? =0 ^) 



\q n v-r n u+p r 
i„ iVa 



,1/2 



in 



,1/2 



1=0 



\1n 



1/2 



Noting that q n € and that g n by Lemma 3.20 completes the proof. 
Corollary 3.22. If h e K \ H(Q), i/ien |g n | tends te oo. 



□ 



Proof. This follows almost immediately from the fact that there are only finitely 
many rational points (J, |) € S that are written lowest terms, are inside the unit 
sphere, and have \q\ bounded. Since the volume of e-radius balls centered at these 
points shrinks to zero as e shrinks to zero, no irrational point h can be arbitrarily 
well approximated by such points. □ 



As a corollary to the proof of Theorem 3.21 we obtain a new form of the classical 
formula for regular continued fractions: 



Pn 



1 



q n (q n +i + q n ■ T n+1 x) ' 



The left-hand side of this formula may be considered to be the distance between x 



and the point p n /q n . Recall that in Theorem 3.21 we showed that 
(3.8) d(K{ 7i }^ ,h) 



Theorem 3.23. Let h € H with continued fraction digits CF(h) — {^i}, asso- 
ciated to a fundamental domain K with vad(K) < 1, and rational approximates 

K{7i}" =1 = ■ Then, in the notation of Lemma 



3.16 



1 

<7^(<7n+l + q n +iu n+ i - q n v n +i) 



Proof. This follows immediately from (3.7) and Lemma 3.19 



□ 



3.5. Uniform convergence. We continue with the assumptions of Theorem 1 3. 2 1| 
and the notation of Lemma |3.16| The purpose of this section is to show that 
the terms {jf^i^^j a ^ so converge to (u,v). There are several ways to prove this, 
including along the lines of the proof of Theorem |3.21[ or using the geometric 
duality between the columns of Q n . The proof below provides some lemmas of 
independent interest. 

Lemma 3.24. We have that 

V 2 k«l < Wn\ < V 2 Wn\\q n -l\ 
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Proof. We associate to h a sequence of unitary transformations A li £ U(2, l;Z[i]). 
We then have Q n = A 7l ■ ■ ■ A ln £ U(2, l;Z[i]). Recall further that for any matrix 
M £ f/(2, 1) we have M^JAO = I. In particular, [7(2, 1) is closed under trans pose 
conjugation, so we have Q^ n £ U(2, l;Z[i]). This matrix has (by Lemma 3.16) the 
form 



q n Pn 
-q n -i —r n -i 




Since Q\ £ U(2, l;Z[i]), we must have that the first column is a null vector with 
respect to J. That is, by (2.3) we have 



\q n \ 2 +2Re(g„g„_i) 







and thus 



,* |2 n / Qn— 1 

■2 5n Re 



Since q n and are non-zero, the real part of q n -\/q n must be non-zero and so 
must have norm at least l/|g„|. Also, the absolute value of the real part of q n -i/q n 
must be at most \q n —i/q n \- This completes the proof of the inequality. □ 



In particular, Lemma 3.24 shows that \q n \ tends to infinity as n grows. 
Lemma 3.25. We have 

'n , qn 7 Pn Pn ^n 



qn qn qn ' qn 



q n q n q n ■ q n 



Proof. Lemma |3 . 1 7| implies that 

q n = -q n r n + f n q n . 

Solving for f n and dividing by q n gives the first equality. 
In addition, Lemma |3 . 1 7| implies that 

Tn PnQn PnQn' 

Solving for p n and dividing by q n gives the second equality. 



□ 



Theor em 3 .26. Under the assumptions of Theorem \3.21\ and with the notation of 
3.16, the points (j^, converge to h — (u,v), provided h is irrational. 



Lemma 



Proof. Since we know by Theorem 3.21 that (r n /q n ,p n /q n ) converges to (u,v), it 
suffices, by Lemma [3.25| to show that 



q n 
q n ■ q n 



and 



Qn ■ Qn 



converge to 0. This is easy to see for the first number, since \q^/q n \ — 1 and \q n \ 
grows to infinity. For the second number, we have that \r^/q n \ — \f n /qn\ which 
converges to and hence is bounded. Since, again \q n \ grows to infinity, the 
second number must converge to as well. □ 



24 



A. LUKYANENKO AND J. VANDEHEY 



4. Heisenberg Gauss Map 

Let K be a fundamental domain for H(Z) with rad(i^) < 1 and T : K — > K the 
associated Gauss map (for this section, we will ignore the fact that T is not well- 
defined at the origin). The purpose of this section is to prove the existence of a 
measure \x on K, absolutely continuous with respect to Lebesgue measure A onR 3 , 
for which T is invariant, and to show that T is ergodic; that is, T~ l E = E implies 
/i(E) — or fj,(K \ E) = 0. We will use the following variant of Renyi's theorem, 



whose terminology we will describe below (see the proof of Lemma 4.9 1 



Theorem 4.1 (See Theorems 4 and 8 in [15]). Let T give rise to a fibred system 
over a set K, with digit set V. Let A be some measure. Suppose 

(1) \{K) = 1; 

(2) All cylinder sets are full; 

(3) For any infinite admissible sequence w — {wi, W2, ■ ■ ■ } of digits from T>, we 
have 

lim diamCi w w w \ =0; 

n— >oo 

(4) There is a constant C > 1 such that for all finite admissible strings w of 
length n, 

su PyeT™c lu JyT n 



mi veTnCm J y T r ' 



< C. 



Then T is ergodic and admits a finite invariant measure [i, which is absolutely 
continuous with respect to A (and vice-versa). 

Remark 4.2. In literature on continued fractions, it is customary to consider 
the inverse transformation to T n given by T~ n (h) = O'fi ■ ■ ■ v^ n h. In the fourth 



condition of Theorem 4.1 one then bounds the Jacobian of T 71 , with notation 

w(7i, • • • ,7«;y) := J y T~ n . 

We will not use this notation since we do not manipulate the digits 7». 

Remark 4.3. The Gauss Map T does not satisfy the conditions of Theorem |4.1| 
Instead, we will define a related map S that does satisfy the conditions. However, 
we start by studying T. 

Let w — {71 , • • • , 7„} be a sequence of elements of H(Z) , and consider the (possibly 
empty) cylinder set C w of points in K whose continued fraction expansion starts 
with w. A cylinder C w is full if T n C w = K, that is, if by removing the first n 
digits from points in C w , we see the full range of points we would expect to see. 
Note that, for example, the cylinder Cm 0)} is not full for any choice of K. For a 
point h £ K with CF(h) = {ji}, let n(h) denote smallest number n such that the 
cylinder C{ 71 ... j7re j is a full cylinder. 

Working in the geometric model of H, let A denote the Lebesgue measure on M 3 = H. 
Note that A is also the Haar measure for H, that is, left multiplication by an element 
of EI does not distort A. 

In the next four lemmas, we study the cylinders C w . We show that as w grows, 
the cylinder becomes smaller and, at some point, it becomes full. Furthermore, the 
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transformation T distorts A on C w in a uniform way for all cylinders at the same 
rank (length of w). 

Lemma 4.4. Let h G K be an irrational point contained on the interior of K and 
CF(h) — Denote the cylinder sets associated to h by C n = C{ 7lj ... l7 „}- 

For sufficiently high n, we have C n C K. Furthermore, 

lim diam(C*„) = lim A(C n ) = 0. 

Proof. By Theorem |3.21| we have that the rational approximates of h converge to 
h. In particular, the rational approximates K{7 i }™ =1 are eventually contained in 
K, and therefore in the associated cylinder set C„. Furthermore, our bounds on the 
Euclidean distance (in the planar Siegel model) from K{7i}™ =1 to any point of C n 
are independent of the particular point. Thus, we have that for sufficiently large n, 
the closures C n of the cylinder sets form a nested sequence of compact sets whose 
intersection is equal to {h}. This implies that the sets C n , in the limit, become 
trivial both in diameter and measure, as desired. □ 



Lemma 4.5. Under the assumptions of Lemma 4-4 < there exists n(h) S N such 
that the cylinder set C n ^ is full. 

Proof. Recall that a cylinder set C n is full if T n C n = K. Recall that for a point 
in Cat, the map T n is given by T n (h') = 7,7 1 t- ■ M . Thus, C n is full exactly 
if C n — L-fi ■ ■ ■ Lj n K. By definition of C„, we have C n = K n L'fi ■ ■ ■ i"/ n K, so it 
suffices to show that for high n the set 171 •• • ij n K does not contain points outside 
of K. 



Now, by Lemma |4.4[ for high n, C n is contained on the interior of K. Furthermore, 
it is easy to see that the set 171 • • • ij n K is connected. Thus, for sufficiently high n, 
471 • • • ij n K cannot have points outside of K, and C n must be a full cylinder. □ 

Recall that the Jacobian Jhf of a measureable transformation / measures the 
volume distortion caused by / at the point h. 

Lemma 4.6. Let h be an irrational interior point of K with CF(h) = {7^°^. Let 
w = {71, . . . ,7n} be a finite subsequence of the continued fraction digits. We then 
have 

n-l 

J h T n = J] \\T(h)f S . 

i=0 

Proof. The transformation is given by T n = r y^ 1 L---^ 1 L. Left translation by 
7i acts by a shear on H in the geometric model, and therefore does not distort 
Lebesgue measure. It suffices to show that for any point h E K , 

J h b = \\h\\~ . 

Recall that by Lemma pO] we have 

d(h,k) 



d(ih, ik) 



\h\\\\k\\ 

Thus, the image of a b all of radius e around h is, in the limit, a ball of radius 



h\\ e. By Lemma 2.3 this implies a distortion factor of \\h\\ . □ 
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Lemma 4.7. Let w = {71, • • • , 7„} be a sequence of elements o/H(Z). Then, 

su PheC J h T n ( 16rad(in 

< exp ' 



mf h £ Cw J h T n ~ ^Vl-rad(if) 2 
Proof. Let h,k G (7 W . By Lemma |4.6| we have 



J h T n /Vr 1 II^A;| 



n 



By the triangle inequality we have: 



\ Tik \\ < 1 , d^h^k) 



lir^H - \\T i h\\ ' 

Furthermore, for i =/= n, 

dCFh^k) _ diiT'h^T'k) _ d(T l+1 h,T l+1 k) 
\\T l h\\ \\tT l k\\ \\lT%\ 



< d(T^KT^k)^d(K) < d[T ^l^ k) rad(X) 2 . 



Continuing inductively, we have that 
d{T i h,T i k) 



\\T% 

Returning to the Jacobian, we have 

j h T n ( n -f} llT^ 1 ^ s " 1 



< d(T n h,T n k)Y&d(K) 2 ( n ~ i '> < 2rad(X) 2( "- l)+1 . 



\2(n-i) + l 



JkTn V n^4 <n(l+2rad(^ 

K \j=0 11 11 / i=0 

/ Ti- 1 \ / n x 

< exp I 8 J] 2rad(/0 2(n_l)+1 < cxp 18^ 2rad(if) 2i+1 

V i=0 / V *=1 / 

' 16rad(X) 



< exp 16 ^ rad( J ftT) 2l+1 J = cxp ( - 



rad(i4T) 2 

as desired. □ 

Definition 4.8. Let the jump transformation S : K — >• K associated to T be 
defined by S(h) = T n ( h \h). Recall that n(h) is the smallest integer so that C n is 
a full cylinder. 

By Lemmas |4.4| and |4.5[ S is well-defined for A-almost-every point of K . 

Lemma 4.9. The jump transformtion S is ergodic and admits an invariant measure 
As, which is absolutely continuous with respect to A. 

Proof. It suffices to show that the jump transformation and its associated collection 
of cylinders satisfy the conditions of Theorem |4.1| 



we 



The map S is defined on the fundamental domain K of H(Z). By Lemma 2.17 
have that \(K) = 1. 

The fibred system associated to 5* arises from the collection of full cylinders C' n ^ 
for points h £ K (these are cylinders with respect to T) . (Note that if k € C n ^ , 
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then nik) = n(h) and so C n ^) = Cn{h)-) O n these cylinders, the maps T n ( h ' and 
S agree, so the cylinders are also full with respect to S. 

By Lemma [4. 4[ the diameters of the cylinder sets C w limit to as the length of w 
increases. 

For the last condition of Theorem |4.1[ we may apply Lemma [477] While the lemma 
is stated for T, it also holds for S since S(h) = T n( * h \h). □ 

Theorem 4.10. Let K be a fundamental domain for H(Z) with rad(i^) < 1, and 
T : K — > K the associated Gauss map. Then T is ergodic and admits an invariant 
measure At, which is absolutely continuous (in fact, equivalent to) Lebesgue measure 
A. 



Proof. In Lemma 4.9 we used a variant of Renyi's Theorem (Theorem 4.1 ) to show 
that the jump transformation associated to T is ergodic and admits an invariant 
measure As which is equivalent to A. 

Let Do = K, and for each n > let D n denote the union of non-full cylinders of 
rank n: 

D n = \jC w , w = {71, • • • , 7„} and C w is not full. 

Theorem 11 in |15) states that T is ergodic with respect to an invariant measure 
At given by 



Ar(£)^A s (r"EnJJ„). 



ra=0 



Since we have that As is equivalent to Lebesgue measure A, it follows directly that 
the same is true for At- □ 



In the case of classical continued fractions, the density of the invariant measure is 
given by a relatively simple expression. It's not clear whether this is the case for 
Heisenberg continued fractions, but it seems likely that the density is, at the very 
least, given by a real-analytic function. One may attempt a proof along the lines 
of Theorem 6.1.3 of [12]. 

As an application of the ergodicity of T, we can attempt to visualize the invariant 
measures on Kc and Kr>. Specifically, the Birkhoff Ergodic Theorem states that for 
a generic point h 6 K, the time average of T n (h) estimates the associated invariant 
measure. The results are displayed in Figure |1J 
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